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1. Introduction 



As so many natural laws and models of natural phenomena are described by 
nonlinear dynamical systems dyi/dt = F(x), it is no exaggeration to say that they 
abound in modern science and technology. Physical and numerical experiments 
show that deterministic chaos in such systems is ubiquitous. On the other hand, 
the recent observations concerning possibility of controlling chaos [1] indicate its 
new promising applications in engineering and medical sciences [2] . No wonder that 
the interest in the theory of chaotic dynamical systems is steadily increasing. 

The object of the present paper is to formulate a hypothesis concerning "quantal" 
Hilbert space criterion of chaos in nonhnear systems dx/dt — F(x), where F is 
analytic in x, with the use of the formalism developed by the author [3-9] , relying 
on reduction of nonlinear dynamical systems to the linear, abstract, Schrodinger-like 
equation in Hilbert space. Namely, it is suggested that the chaotic behaviour of the 
nonlinear dynamical systems can be related to the growth of the "quantal" entropy 
which can be naturally introduced within the Hilbert space approach resembling 
quantum mechanics. 



2. Hilbert space approach 



In this section we briefly outline the Hilbert space description of nonlinear ordi- 
nary differential systems [3] . Consider the analytic system (complex or real) 

d7 

- = F(z), z(0) = Zo, (1) 

where F: C'^ C'^ is analytic in z. 
The vectors \z,t) defined by 

|z,t)=exp[|(|z|2-|zoP)]|z), (2) 

where |z) is a normalized coherent state (see appendix) and z fulfils (1), satisfy the 
following linear, Schrodinger-like equation in Hilbert space: 

j^\z,t)^M\z,t), 1^,0) = |zo), (3) 
where M is the boson operator such that 

M = a^.F(a). (4) 
In view of (2) we arrive at the following eigenvalue equation: 

a|zo,t) = z(zo,t)|zo,t) (5) 
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relating the solution z(zo,t) of (1) and the solution |zo,t) of (3). Thus, it turns 
out that the nonlinear dynamical system (1) can be cast into the linear, abstract, 
Schrodinger-like equation (3). We note that the algorithm can be immediately 
extended to the case of nonautonomous systems such that the corresponding vector 
field is analytic in z-variables. We also remark that the eigenvalue equation (5) 
suggests the "quantization scheme" of the form z — > a, where z fulfils (1), within 
introduced formalism resembling quantum mechanics. 

Observe that (3) and (5) correspond to the "Schrodinger picture" within the 
actual "quantal" treatment. We end this brief account of the Hilbert space approach 
with discussion of the "Heisenberg picture" . The "Heisenberg equations of motion" 
obeyed by the time-dependent Bose annihilation operators are of the form 

^ = [a,M], a(0)=a, (6) 

where M is the "Hamiltonian" given by (4). Since the "Hamiltonian" in "Schrodinger 
picture" M is independent of time, therefore it coincides with the "Hamiltonian" in 
"Heisenberg picture" . 

The formal solution of (6) can be written as 

Si{t)^V{t)-'aV{t), (7) 

where V{t) — e*^ is the evolution operator. It is clear that the following relation 
holds true: 

\zo,t) ^V{t)\zo), (8) 

where |zo,t) is the solution of (3). 

Now eqs. (5), (8) and (7) taken together yield 

a(t)|zo) = z(zo,t)|zo). (9) 

Hence we find that the solutions of the systems (1) coincide with the expectation 
values (covariant symbols) of the time-dependent Bose annihilation operators, i.e. 

z(zo,t) = (zo|a(t)|zo). (10) 

It should be noted that whenever the "quantization scheme" z — a is assumed, 
where z fulfils (1), then (10) forms the "Ehrenfest's theorem" within the actual 
"quantal" approach. 



3. The criterion of chaos 



We now formulate the hypothesis concerning the criterion for the nonlinear dy- 
namical systems (1) to show chaotic behaviour. Consider the system (1). We intro- 
duce the operator p of the form (see (A.2)): 

p = y dAi(w)e-l*-^°lV)(w|. (11) 
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Evidently, the operator p is Hermitian and nonnegative, i.e. 

= p, (12a) 
p > 0. (12b) 

Furthermore, using (A. 6) we immediately obtain 

Trp = 1. (12c) 
Finally, taking into account (11), (A. 2), (A. 5) and (10) we get 

(a(t)) = Tr(pa(t))=z(t), (12d) 

where a(t) are the time-dependent Bose annihilation operators and 2,{t) fulfils (1). 

We have thus shown that if the "quantization scheme" z — a, where z satisfies 
(1) is assumed, then p plays the role of the density matrix within the "quantal" 
Hilbert space approach. Prompted by this analogy, we define the "entropy" as 

= -Tr(plnp). (13) 

On using the relation 

Trp" = (2^7^' ^ ^ 1' (14) 
following directly from (11), (A. 2) and (A. 4), we find 

5 = -vf' + ^~l^-i-ln-l- = 2;tln2. (15) 

Since an arbitrary vector in a Fock space can be specified by an infinite sequence 
of n- vectors of I's and O's and the number of different n- vectors is 2^*, therefore 
the obtained value of the "entropy" can be regarded as an averaged amount of the 
information necessary to fix the coherent state |zo) in a Hilbert space of states. 

Now in analogy to quantum mechanics, we introduce the time-dependent "den- 
sity matrix" such that 

p{t)^ j dp(w)e-l"'-"°l'e*^|w)(w|e*^\ (16) 

where M is the"Hamiltonian" related to the system (1). 

We note that p{t) is Hermitian and nonnegative at any time. The time-dependent 
"entropy" corresponding to (16) is given by 

S{t) = -T,[p{t)\np{t)]. (17) 

We are now in a position to propound our hypothesis. We assert that if the following 
conditions hold: 

I V |z(zo,i)r<C, (18a) 
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where z(zo, t) is the solution to (1); C > is a constant, 

V ^ > 0, (18b) 
t>u dt ^ ^ 

then the system (1) is chaotic one. 



4. Discussion 

We have formulated in this work a hypothesis concerning a new, "quantal" 
Hilbcrt space criterion of chaos in nonlinear dynamical systems. Note that the 
condition (18a) ensures that the "density matrix" p{t) is of trace class for all t > t*. 
Indeed, from (16) and the relation 

{zo,t\zo,t) =exp(|z(zo,f)|^ - |zo|^), (19) 

which is an immediate consequence of (2), it follows that 

Trp(i) = y (i//(w)exp(-|w-zo|^)exp(|z(w,t)|^ - |w|^), (20) 

where z(zo,t) is the solution to (1). 

On the other hand, the condition (18a) is consistent with the exceptional role 
of dissipativity of systems showing chaotic behaviour. We remark that the nonuni- 
tary evolution implied by the Schrodinger-like equation (3) is crucial for the ac- 
tual treatment. In fact, as in quantum mechanics, the unitary evolution implies 
S{t) = S = const. We note at the same time that the evolution operator V{t) = e^^ 
defined by (8) is unitary only in the case of the linear system (1) with a skew- 
Hcrmitian fundamental matrix. Finally, we point out that the "quantal" Hilbert 
space approach recognizes the divergence of orbits in the phase space of the system 
(1). Namely, we have the formula 

|(zo,t|Zo + (5zo,t)P , . , , r .N f /o1^ 

= exp(-|z(zo + ()zo,i) -z(zo,i)r) (21) 



(zo, t|zo, t) (zo + 5zo, t\zo + 6zo, t) 



which can be easily obtained from (2) and (A.l). Therefore, the growth of the 
distance between the two trajectories of the system (1) starting from zq and zo + 5zo 
implies the decay of the corresponding transition probability from the left-hand side 
of (21). 



Appendix. Coherent states 
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We append the basic facts about coherent states. Recall first that the Bose cre- 
ation (a''^) and annihilation (a) operators, where a^ = {a\, . . . , a|.), a = (ai, . . . , a^), 
obey the Heisenberg algebra 

[tti, a]] = SijI, 

[ai, Qj] = [a], a]] = 0, i, j ^ 1, . . . , k. 

The coherent states |z), where z e C^, are usually defined as eigenvectors of the 
annihilation operators, that is 

a|z) = z|z). 

The normalized coherent states can be defined as 

|z) = exp(-||z|^)exp(z-a^)|0). 
where |zp = Yli=i l^iP, u-v = J2i=i UiVi and |0) is the vacuum vector satisfying 

a|0) = 0. 

The coherent states are not orthogonal, namely 

(z|w) = exp[-|(|zp + |w|2 - 2z*.w)], 

where the asterisk designates the complex conjugation and z* = {zl, . . . , z^). 
Hence we find 

|(z|w)|2 = exp(-|z-w|2). (A.l) 

These states form the complete (overcomplete) set. The resolution of the identity 
for the coherent states is given by 

rf/x(z) |z)(z| = /, (A.2) 



k ^ 



where 



dii{z) = Y[ —d(ReZi) d{lmzi). 



Now let 10) be an arbitrary state. It can be easily shown that the function (symbol 
of the vector) 0(z*) = (z|0) is of the form 

0(z*) = 0(z*)exp(-||z|2), (A.3) 

where 0(z*) is an analytic (entire) function. 
Taking into account (A.2) and (A.3) we get 

J rf/^(z)exp(-|zn(0(z*))>(z*). 

R2fe 
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Thus the abstract vectors can be represented by analytic functions. This represen- 
tation is usually known as the Bargmann representation. The action of the Bose 
operators in the Bargmann representation has the following form: 

a0(z*) = at0(z*) = z>(z*). 

On using (A. 2) and (A. 3) we arrive at the following reproducing property of coherent 
states: 

0(w*)= y dA*(z)exp(-|z|2)3C(w*,z)0(z*), (A.4) 

where the reproducing kernel (Bergman reproducing kernel) is 

3C(w*,z) = exp(w*'z). 

Taking the Hermitian conjugate of (A.4) we obtain the following form of the repro- 
ducing property: 

^{w) = J (i/i(z) exp(— |zp) exp(w-z*)-?/'(z). (A. 5) 

We finally remark that the coherent states are the convenient tool for the study of 
operators. For example, the trace of a linear operator L can be expressed by 

T^L^ J (i//(z)L(z*,z), (A.6) 
where 

L(z*,z) = (z|L|z) 
is called the covariant symbol of the operator L [10]. 
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